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1. Introduction. The diffuse component of the galactic high-energy (£50 
MeV) gamma rays results mainly from the interaction of CR nuclei (mostly 
protons with energies of a few GeV) and electrons (with energies up to 
several hundreds of MeV) with the nuclei of the interstellar gas (via the 
decay of tt "- mesons and bremsstrahlung) . An additional contribution is 
obtained from the interaction of CR electrons (with energies > 1 GeV) with 
the interstellar photons (mainly in the optical and infrared range) 
through the inverse-Compton (IC) process. Gamma-ray astronomy therefore 
offers an excellent means to study the distribution of CR particles 
throughout the Galaxy, but it is essential to know the distribution of 
the target interstellar gas particles, the major constituents being 
atomic (HI) and molecular (H 2 ) hydrogen. Although large-scale mapping of 
the HI component had been performed using its characteristic 21-cm line, 
the numerous previous gamma-ray studies of the galactic CR distribution, 
using gamma-ray observations obtained by the SAS-2 and COS-B satellites, 
suffered severely from uncertainties in the galactic distribution of 
interstellar molecular hydrogen. Large-scale millimetre-wave surveys of 
the CO molecule covering more than half of the Milky Way, obtained with 
the Columbia 1.2m telescopes, are currently available and can be used to 
trace the H 2 ; the COS-B observations have sufficient resolution and 
sensitivity to constrain the relation between the integrated CO line 
intensity W qq and the molecular-hydrogen column density N(H 2 ) (Lebrun 
et al., 1983; Bloemen et al., 1984a). 

The velocity information of the HI and CO observations is used as a 
distance indicator to ascertain the spatial distribution of the 
interstellar gas. Using this distance information, the galacto-centic 
distribution of the gamma-ray emissivity (the production rate per H atom) 
is determined for three gamma-ray energy ranges (70 MeV - 150 MeV, 150 
MeV - 300 MeV, and 300 MeV - 5 GeV) from a correlation study of the 
gamma-ray intensity maps and the gas-tracer maps for selected galacto- 
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centric distance intervals, taking into account the expected IC contribu- 
tion and pointlike gamma-ray sources (details are given in Section 2). On 
the assumption that unresolved gamma-ray point sources do not contribute 
significantly to the observed gamma-ray emission, the gamma-ray 
emissivity is proportional to the CR density and, more specifically, the 
energy dependence can be used to study separately the distribution of CR 
electrons and nuclei: whereas the emission for the 300 MeV - 5 GeV range 
is dominated by TT°-decay, the 70 MeV - 150 MeV range has a large electron 
bremsstrahlung contribution. For a complete description of the method and 
the results the reader is referred to Bloemen et al. (1985). 

Bloemen et al. ( 1984b, c) applied such a method to the C0S-B 
gamma-ray observations and HI 21-cm line observations of the second and 
third galactic quadrants alone (they showed that H 2 can be neglected in 
the outer Galaxy within the uncertainties of the analysis), and found 
evidence that the distribution of the CR electron density decreases 
beyond the solar circle while the density of the CR nuclei is approxi- 
mately constant out to large (-20 kpc) galacto-centric distances. 

2. Method. Skymaps of HI column densities, N(Hl)i, and integrated CO 
line intensities, W qq £, were constructed in four (i - 1,2, 3, 4) 
galacto-centric distance ranges (2 kpc < R <8 kpc, 8 kpc < R < 10 kpc, 10 
kpc < R <15 kpc, and R > 15 kpc) and convolved with the energy-dependent 
COS-B point-spread function. We selected these four distance intervals 
because the angular distributions of the gas in each interval show 
distinct differences, which are needed to ascertain the contribution of 
the gas in each interval to the observed gamma-ray intensities. 
Similarly, determining the conversion factor between Wqq an< * N(H2) 
(i.e. the ratio X = N(H 2 )/W C q) requires distinct differences between 
the structures in the HI and CO maps. Due to the limited sky area that is 
covered by the CO observations, the following correlation analysis has 
been restricted to the first and second galactic quadrants and the Carina 
region (270° < 1 < 300 0 ) , and the latitude range -4. 5° < b < 6 . 5 0 . 

We investigated which combination of gamma-ray emissivities , q£j 
(j = 1,2,3 corresponds to the 70 MeV - 150 MeV, 150 MeV - 300 MeV, ana 
300 MeV - 5 GeV ranges), best describes the observed gamma-ray 
distributions, Iv assuming a relation of the form: 

4 ,J 

I Y ,j "{ £ iy_-(N(Hl)i + 2Y j .W c0 , i ) )+ lie, j + *b f j- 
* i= 1 4 tt 

The term enclosed by braces represents the intensities that 
originate from CR collisions with atomic and molecular hydrogen, the term 
IlC,j represents the modelled (small) IC contribution (Bloemen, OG3.1-2), 
and* 1^ j is the total isotropic gamma-ray background, including the 
(dominant) instrumental background. The parameters Yj are equal to the 
ratio X, independent of energy, if CR particles are not excluded from, or 
concentrated in, molecular clouds. Pollock et al. (1985) have shown that 
some gamma-ray point sources are present in the first quadrant; these 
have been included in our model as described by Pollock et al. 

We applied a likelihood analysis on l°xl° bins to determine the 
values and formal statistical uncertainties of the parameters q^ j , as 
well as of Yj and Ib,j* Starting from the general model with these six 
parameters for each ’energy range, we tested (using the likelihood- 
ratio) whether various simpler models with fewer parameters (i.e. 
constant emissivity distributions as a function of R, identical Yj values 
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Figure 1: Galacto- centric distribution 
of the gamma-ray emissivity for three 
energy intervals j. The crosses s 
together with statistical 1 a error 
bars j indicate the fit values for 
selected distance ranges i ( indicated 
by the vertical lines ) without any 
constraints. The circles are the 
emissivity values when the values of 
Yj are forced to be identical. The 
full curves represent the emissivity 
distributions Qj(R) in the case of 
exponential distributions for the CR 
electrons and nuclei. The dashed 
curves indicate the ^ °-decay contribu- 
tion from the CR nuclei (f = 1). 


and/or identical emissivity distribu- 
tions for each energy range) give 
significantly worse fits to the data. 


3. Results and discussion. Figure 1 shows the resultant radial 
emissivity distributions. Since neither the two emissivity values outside 
the solar circle, nor the two values inside, can be determined entirely 
independently, these radial distributions have to be judged carefully. 
The likelihood-ratio hypothesis testings showed: (1) a gamma-ray 
emissivity gradient as a function of galacto— centric radius is required 
for the integral 70 MeV — 5 GeV range, (2) the gamma-ray model displays 
an energy dependence, (3) this energy dependence can be accounted for 
fully by an energy-dependent emissivity gradient (so energy- independent 
Yj values) that is strongest for low energies, (4) it can less satisfac- 
torily be ascribed to different emissivity spectra of HI and H 2 (i.e. 
different Yj values). The ratio X = N(H 2 )/W^q was found to be 
(2.75 ± 0.35)xl0~ 26 mol. cm -2 K -1 km -1 s, and additional tests showed 
that X is constant throughout the Galaxy, within uncertainties. This X 
value is independent of excitation and abundance effects, which have 
plagued previous determinations. It should strictly be regarded, however, 
as an upper limit if a population of unresolved galactic gamma-ray 
sources distributed like CO exists. The resultant H 2 mass is equal to the 
HI mass for 2 kpc < R <10 kpc, namely 0.9x10^ M Q . 

Assuming exponential distributions e s ( R-R 0 ) (R 0 = 10 kpc) for the 
CR electrons (e) and nuclei (n) , the total gamma-ray emissivity at 
galacto-centric radius R can be written as: 


q j (R) = e Se(R R ° ) .{qj(R 0 )-f.q 1I <> > j } + e Sn(R Ro) f. 


qir° 


where the parameters the local tt ““decay gamma“*ray emissivi— 
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Figure 2 : Radial CR-density distributions of the form efHR-Ro) s for 
CR electrons (e) and nuclei (n). The hatched areas encompass 0.8 ^ f^ 1.2 
and 1 o uncertainties for S n and S ej that is - 0 . 21 < S e < - 0 . 09 kpc~ 1 and 
-0.055 < S n < + 0.015 kpe~l. Rote that urithin these areas S e and S n are not 
independent: a stronger electron gradient implies a stronger gradient for 
the nuclei (see table > presenting S e and S n for different f values). 

ties (we estimated these values from the work of Stephens and Badhwar 
(1981), based on the demodulated proton spectrum near the earth) and f is 
an energy- independent scaling factor, because the demodulated proton flux 
may not be typical for the local interstellar medium. Considering a range 
of possible f values (0.8 £f £1.2), we applied a likelihood analysis to 
determine the values of the local emissivity values qj(R 0 ), and of S e and 
S n . Note that we did not fit the emissivity distributions shown in 
Figure 1; our method allows for the dependencies mentioned above (for 
details see Bloemen et al., 1985). The resultant values of S e and S n are 
given in Figure 2 and the corresponding emissivity distributions are 
included in Figure 1. The total local gamma-ray emissivities qj(R 0 ) were 
found to be independent of the f value chosen and are consistent with the 
improved local emissivities determined at medium latitudes (Strong et 
al., OG3.1-3); the implications for the local electron spectrum are 
discussed by Strong (OG3.1-7). The CR electron gradient is required, but 
the results are consistent with a constant density of CR nuclei on a 
galactic scale. The quoted CR gradients could be affected by a population 
of unresolved galactic gamma-ray sources. 

The gradient in the distribution of the CR electrons confirms their 
galactic origin. The exponential scale length of the CR nuclei distribu- 
tion from the present analysis is at least 15 kpc; this is much larger 
than for the type of objects generally considered as candidates for CR 
sources (e.g. supernovae and pulsars, both having a scalelength of ^5 kpc 
(Kodaira, 1974; Lyne et al., 1985)). The gradient of the CR nuclei, if it 
exists, i9 so weak, that on the basis of gamma-ray observations, it can 
no longer exclusively be claimed that CR nuclei (with energies of several 
GeV) are produced in the Galaxy. If the latter is the case, and if they 
are produced by objects like pulsars and supernova remnants, then the CR 
distribution does not seem to reflect the distribution of the CR sources. 
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